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1. Introduction 

Recent observational data indicate that our universe is accelerating. This acceleration 
is explained in terms of the so-called dark energy (DE), which may be explained in 
modified gravity models (for a general review see Ref. pQ). DE could also result from 
a cosmological constant, from an ideal fluid with a different form of equation of state 
and negative pressure, a scalar field with quintessence-like or phantom-like behavior 
(see [2] and references therein), etc. The choice of possibilities reflects the indisputable 
fact that the true nature and origin of DE has not been convincingly explained yet. It 
is not even clear what type of DE is more seemingly to explain the current epoch of 
the universe. Observational data point towards some kind of DE with an equation of 
state (EoS) parameter which is very close to -1, maybe even less than -1 (the so-called 
phantom case). A quite appealing possibility is the already mentioned modification of 
General Relativity (GR). Modifications of the Hilbert-Einstein action by introducing 
different functions of the Ricci scalar have been systematically explored, the so-called 
F(R) gravity models, which reconstruction has been developed in Refs. [3]- [5]. As is 
known, F(R) gravity can be written in terms of a scalar fieldquintessence or phantom 
likeby redefining the function F(R) with the use of a scalar field, and then performing 
a conformal transformation. It has been shown that, in general, for any given F(R) the 
corresponding scalar-tensor theory can, in principle, be obtained, although the solution 
is going to be very different from one case to another. Also, attention has been paid to 
the reconstruction of F(R) gravity from a given scalar-tensor theory. It is known, too, 
that the phantom case in scalar-tensor theory does not exist, in general, when starting 
from F(R) gravity. In fact, the conformal transformation becomes complex when the 
phantom barrier is crossed, and therefore the resulting F(R) function becomes complex. 
These situations where addressed in [8] in detail, where to avoid this hindrance, a dark 
fluid was used in order to produce the phantom behavior in such a way that the F(R) 
function reconstructed from the scalar-tensor theory continues to be real. 

On the other hand, it has also been suggested in the literature (see Refs. [9], [TTj ) 
to consider modified Gauss-Bonnet gravity, that is, a function of the GB invariant. 
Different cosmological properties of modified gravity models of this kind have been 
studied in Refs. [9]-[25]. Both possibilities have, in principle, the capability to explain 
the accelerated expansion of the Universe and even the primordial inflationary phase 
(see [3] and [8]), with no need to introduce a new form of energy. In this paper we 
will study some specific modified Gauss-Bonnet theories and, by using a technique 
developed in Ref. [5], the corresponding cosmological theory will be reconstructed for 
several cosmological solutions. 

We will here work with the spatially-flat FRW universe metric 

3 

ds 2 = -dt 2 + a(t) 2 Y,(dx 1 ) 2 , (1) 
i=i 

where a(t) is the scale factor at cosmological time t. In GR, the corresponding system 
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of equations are the usual Friedmann equations, namely 

^H 2 = p, l(2if + 3iJ 2 ) = -p, p = -3H(p + p). (2) 

The last of them is the continuity equation for a perfect fluid. As usually, the Hubble 
rate H is defined by H = a/a. In ([2]), p and p are the matter energy-density and 
pressure. The Gauss-Bonnet invariant is given by 

G = R 2 - AR^Ri" + R^xaR^ Xa (3) 

The cosmological models coming from the different versions of modified GB gravity here 
considered will be carefully investigated with the help of several particular examples 
where calculations can be carried out explicitly. 

The paper is organized as follows. In Sect. II, the general technique which will be 
used to reconstruct a gravity theory from a given cosmological evolution setup will be 
explained in detail. Several specific examples will be worked through, for the simple 
case where one considers a Hilbert-Einstein action plus a function on the GB invariant. 
In Sect. Ill, a general model whose action depends on a function of the Ricci scalar and 
of the GB invariant will be studied. With the help of several explicit examples we will 
show how the reconstruction of the cosmological solutions is carried out. In Sect. IV the 
more involved case of a scalar-tensor theory where the GB invariant is coupled to the 
scalar field will be considered. It will be shown there that any cosmological solution can 
be reproduced by an specific scalar potential and a convenient coupling term. Finally, 
in the conclusions we will provide a summary of the main results obtained in the paper. 

2. The [R + /(G)]-model 



S = J d 4 Xy/—g 



i R + f(G) + L r 



(4) 



Consider the following action, which describes General Relativity plus a function of the 
Gauss-Bonnet term (see Refs. [9] and [TT]). 

" 1 

2k 

where k 2 = 87tGjv, Gn being the Newton constant. By varying the action over g pu , the 
following field equations are obtained 

= ^(-^ + l^R) + TV* + \gVf{G) - 2f G RR^ + Af G R^ - 2f G R^R^ T 

-4f G R^R pa + 2(VV7 G )fi - 2g^(V 2 f G )R - 4(V P V^ G )^ - 4(V p V7g)F 
+4(V 2 /g)£T + 4g^(V p VJ G )R pcT - A(V P VJ G )R^. (5) 

For the metric (JT]), these equations give the first FRW equation which has the form 

= -^H 2 + Gf G -f- 2AGH 3 f GG + p m . (6) 
The Hubble rate H is here defined by H = a/a, while G and R are given by 

G = 24(HH 2 + H 4 ), R = 6{H + 2H 2 ). (7) 
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The matter energy density p m satisfies the standard continuity equation 

p m + 3H(1 + w) Pm = 0. (8) 

Let us now rewrite Eq. ([6]) by using a new variable N = In = — ln(l + z), that is, the 
number of e-foldings, instead of the cosmological time t, where z is the redshift (this 
method has been implemented in Ref. [5] for f(R) gravity). The following expressions 
are then easily obtained 

Eq. (jHJ) can thus be expressed as follows 

= -p# 2 + 2AH 3 (H' + H)f G - f 

-576H 6 (HH" + 3H' 2 + AHH') f GG + p m , (10) 
where G and R are now 

G = 2A(H 3 H' + H A ), G = 2A(H A H" + 3H 3 H' 2 + AH 4 H'), 

R = 6(HH' + 2H 2 ). (11) 
By introducing a new function g as g = H 2 , we have 

H = Jg, H>=y^g>, H" = -\g-V 2 g' 2 + l -g^ 2 g". (12) 
Hence, Eq. ( flOj) takes the form 



= -^g+vg(g'+2g)fG-f-24 2 g 



2 // i /2 , <-, 2 / 

2 9 9 + 2" 9 9 



fGG+Pm , (13) 



where we have used the expressions 

G = 12gg> + 24g 2 , G = \2g^ 2 \g 2 g" + gg' 2 + 4^1 , 

R = 3g' + 12g. (14) 
Finally, we can write the FRW equations in a slightly different form [12] 



j^H 2 = p eff , 1(2H + 3H 2 ) = - Peff , (15) 
where the effective energy and pressure densities are 

Peff = WeffPeff, Peff = -3H(p eff + p eff ), (16) 



with 



Here 



Peff = pG + Pm, Peff =PG+Pm- (17) 



PG = Gf G -f- 2AH 3 Gf GG , 
Pg = -Pg + 8H 2 G 2 f GGG 

-192f GG (AH 6 H - 8H 3 HH - 6H 2 H 3 — H A H -3H 5 H - 18H A H 2 ). (18) 

As it was pointed out in Ref. [12], for de Sitter space p G = w G p G = —p G , that is 
w G = -1. 
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2.1. Example 1 

As a first example we consider the ACDM-model. As will be shown below, for the 
gravity theory described by the action (J4j), it is possible to reconstruct such evolution 
with no need of a cosmological constant term. For the ACDM case, we have 

g = H 2 = Hi + *^V» = H 2 + ^a^e~™ = H 2 + /*, (19) 

where I = k p ° a ° , z = e~ m . Hence, we get 

z = -3z, g = -3lz, g" = 9lz. (20) 
Finally, for G we obtain the expressions 

G = 24H*+12H 2 lz-12l 2 z 2 , G = -3zHW = -3Hz(12H 2 l-24l 2 'z), (21) 

where W = 12H^l — 2Al 2 z. We can reverse Eq. (12T1) and the number of e-foldings can 
be written as a function of the Gauss-Bonnet term. This yields 

6/ 3 3H 2 ± ^81# 4 - 3G 

Hence, it follows that 



_ 3H 2 ± - 3G 9H 2 ±J81H$-3G 
g = H 2 = H 2 + ° V 6 = ° V 6 . (23) 

Finally, for the function f(z) we get the following equation 

a 2 f zz + aj z + a f + b = 0, (24) 

where 

a = -w 2 = -24 2 /V + 24 2 H 2 l 3 z - 144ff 4 / 2 , 

oi = 12(H 2 + lz)[U4l 2 z(H 2 + Iz) + (2H 2 - lz)w], 

a 2 = 72zw(H 2 + lz) 2 , b = [p m - ^(H 2 + lz)]w 2 . (25) 
The energy density can be expressed as 

p m = -Dz - E, (26) 

being 

D = UAHHtl^-- + —l E = UAHHtl 2 [--^- 2 --]. (27) 
[ 5k 2 9 9 J [ 9 5k 2 9 J V ; 

Here 5 = const. As a consequence, Eq. (|24|) has the following particular solution 

f{z) = dz 2 + $z + H 2 5, (28) 

where 

6, = ^" a -^ + 1 )]' ^ = 4*" 2 + ^ + 1)] - (29) 
Thus, we have 



3H 2 ± J81# 4 - 3G , 3# 2 ± J8m„ 4 - 3G 
/(G) = fl[ ° V 6/ ] 2 + 0[ ° V ^ ] + H 2 6. (30) 

We observe that this function reproduces exactly the same behavior as the ACDM model 
in the context of Gauss-Bonnet gravity. 
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2.2. Example 2 

Let us now consider a second example 

g = H 2 = le- 3N = Iz. (31) 
This case describes a time evolution given by 

H(t) = -±- , (32) 
t — 1 

which is equivalent to the cosmological evolution of a pressureless fluid in GR, what 
gives a decelerated expansion. As in the example above, we find that z = ±i^p-. In 
the present case, we have the following differential equation, 

*7« - \zf z + \f - C - = 0, (33) 
6 3 3 

where c = p m - ^lz. 

Let us assume that c = 0, then Eq. (1331) takes the form 

z 2 f zz - -zf z + \f = 0. (34) 
o 3 

This equation admits the following exact solution 

h{z) = C 1 z 2 + C 2 z*. (35) 
As a consequence, for this second example we obtain the following model 

f(G) = -Ci^ + C 2 (^)iG^ = C[G + C 2 G&. (36) 

It is now easy to see from this example that the GB terms could actually contribute in 
the matter dominated epoch of the universe evolution, provided the f{G) function is 
taken properly, what is certainly an interesting result. 



2.3. Example 3 

As third and last example we will now consider the case when the Hubble parameter 
and the matter energy density behave as 

tt _ mN n _ i , 3 2mN , 96(m + 1) 5mN 

tl — e , p m — bi + — e H o 2 e , 

where bj = const. This case corresponds to the so-called phantom behavior, what means 
that we have an effective EoS parameter w e ff < — 1, and produces a superaccelerated 
phase that could end in some kind of future singularity (see Refs. [25] and [26]). We 
first write the Hubble parameter as a function of time, 

H(t) = ^ , (38) 

where t s , usually called the Rip time, defines the moment when the Big Rip singularity 
occurs. We can find the f{G) function that reproduces this model by solving the 
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equation ( TTU1) . This yields following solution 

F(G) = F(N) = h + 96(m 5 + \ e 5mjV = 
96(m + l), r G .5 

3. The F(i?, G) model 

Let us now consider a more general model for a kind of modified Gauss-Bonnet gravity. 
This can be described by the following action 

1 



S = J d A x^[^F{R, G) + L m \. (40) 
Varying over g^ u the gravity field equations are obtained [12] 



q _ rppv _|_ _ 2F G RR flu + 4FcRpR up 

72 



-2F G R waT R v paT - \F G R wav R pa + 2(V P V U F G )R - 2g^(V 2 F G )R 

-4(V p V^F G )R vp - 4(V p V u F G )Rw + 4(V 2 F G )i?^ + 4^ iy (V p V (J F G ) J R^ 

-4(V p V (J F G ) J R wl/,J - F G RT + V^V'Ffi - g^V 2 F R . (41) 

In the case of a flat FRW Universe, described by the metric (jTJ), the first FRW equation 
yields 

= ^(GF G -F- 2AH 3 F Gt ) + 3(H + fT 2 )F fl - 3ifF ra + k 2 p m . (42) 

And from here, using the techniques developed in the above section, it is plain that 
explicit F(R, G) functions can be reconstructed for given cosmological solutions. 

3.1. De Sitter Solutions 

As well known, the de Sitter solution is one of the most important cosmological solutions 
nowadays, since the current epoch has been observed to have an expansion that behaves 
approximately as de Sitter. This solution is described by an exponential expansion of 
the scale factor, which gives a constant Hubble parameter H(t) = H . By inserting it in 
the Friedmann equation ()42l) . one finds that any kind of F(R, G) function can possibly 
admit de Sitter solutions, with the proviso that the following algebraic equation has 
positive roots for H 

= \(G F G (G ) - F{G Q , Rq)) + 3H 2 F R (R ) , (43) 

being Rq = 12Hq and G = 2AHq] we have here neglected the contribution of matter for 
simplicity. As it was pointed in Refs. [7] and [S] for the case of modified F(R) gravity, 
the de Sitter points are critical points for the Friedmann equations, what could explain 
the current acceleration phase as well as the inflationary epoch. This explanation can 
be extended to the action (HOI , so that any kind of function F(R,G) with positive real 
roots for the equation (H3j) could in fact explain the acceleration epochs of the Universe 
in exactly the same way a cosmological constant does. 
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3.2. Phantom dark energy 

Let us now explore the cosmic evolution described by ( |37l) in the context of the action 
( |40l) . This solution reproduces a phantom behavior, i.e. a superaccelerated expansion 
that, as pointed by recent observations, our Universe could be in — or either close to 
cross the phantom barrier. We can now proceed with the reconstruction method, as 
explicitly shown in the section above, and a F(R, G) function for the Hubble parameter 
( 1371) will be obtained. For simplicity, we consider the following subfamily of functions 

F(R,G) = f 1 (G)+f 2 (R) . (44) 

Correspondingly, the Friedmann equation ( 1421) can be split into two equations, as 

= -24H*Gf 1GG + Gfic - fi , 

= -3HRf 2RR + 3(H + H 2 )f 2R - l -f 2 + K 2 p m (45) 

For the example (137|) . that is H = e mN , the Ricci scalar and the Gauss-Bonnet 
terms take the following form, 

G = 24(m + l)e imN = 24(m + l)H\ 

R = 6(m + 2)e 2mN = 6(m + 2)H 2 . (46) 
Hence, the first equation in ( 1451) can be written in terms of G, and this yields 

GVlGG _^ G/lG + ^0. (47) 

This is an Euler equation, easy to solve, and yields 

f 1 (G) = C l G 1+2 ^+C 2} (48) 

where Cx t2 are integration constants. In the same way, for the case being considered 
here, the second equation in (j4"5]) . for R, takes the form 

R J2RR ^ Rj2R + ~T. J2 Pm - U . (49) 

2m 2m m 

In absence of matter (p m = 0) this is also an Euler equation, whose solution yields 
f 2 (R) = kiR»+ + k 2 R»- 

l _|_ m+1 4- , /l I ( m + 1 ) 2 m+3 

where /x± = ^ v 4mJ UL- , (50) 

and fci 5 2 are integration constants. Then, the complete function F(R,G), given in fT44|) . 
is reconstructed (in absence of matter) yielding the solutions fT4H|) and fl50|) . The theory 
(3.12) belongs to the class of models with positive and negative powers of the curvature 
introduced in [27]. 

Let us now consider the case where matter is included. From the energy 
conservation equation ([8]) we have that, for a perfect fluid with constant EoS, p m = 
uJmPm, the solution is given by p m = p e~ 3 ^ 1+u ' m ^ iV . By inserting the expression for R 
fl4"6"|) into this solution, we get 

_ 3(l+m m ) 

R \ 2m 
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In such case the general solution for f 2 is given by 

f 2 (R) = hR^+ + k 2 R p ~ + kR A , where, 



2 p m0 (6(m + l)) A 3(1 + w r , 



A(A-l)-l/2m 2m v ; 

Hence, we see that the solution for the Hubble parameter (1371) can be easily recovered 
in the context of modified Gauss-Bonnet gravity. Nevertheless, it seems clear that, 
for more complex examples, one may not be able to solve the corresponding equations 
analytically and it could require numerical analysis. 

4. Reconstruction of scalar-GB gravity 

In this section we consider a four- dimensional action containing the Einstein-Hilbert 
part, a massless scalar field, and the Gauss-Bonnet term coupled to the scalar field. 
The corresponding action is 

S = J d 4 x^g[^R - \d^d^ - V(4>) - £{<f>)G + L m }. (53) 



= — (— + -g^R) + T pv + -d p (f)d v <p - 7 g^d p 
+7;9^(-V + £G) - 2£RBr - 4^R p R up - 2£R ppaT R u p(7T 



An action of this kind has been proposed as a model for dark energy in Ref. [10] . It can 
actually be related with modified GB gravity, as explained in Refs. [13] and [28]. The 
variation of this action over the metric g pv gives [4] 

±(-R^ + l -g^R) + + \d p <\>d^ - 1, 

1 

— ( 

T 

+^R ppU(7 R pcT + 2{V p V u i)R - 2g pv {V 2 £)R - 4(V p V^)R vp 

-4(V P V"0^ P + 4(V 2 0^ M " + Ag pv {V p V (T C)R pa + 4,{V P V a C)R ppvrj . (54) 

In the FRW universe case, these equations and the equation of motion for the scalar 
become [13] 

= -^H 2 + + V{(j>) + 2AH*i + p m , (55) 

= ^(2H + 3H 2 ) + - V{4>) - 8H 2 £ - 16H(H + H 2 )Z + p m , (56) 

= j ) + 3H<P + V; + ^G. (57) 

Note that, as a consequence of the Bianchi identity, Eq. (j57p is satisfied automatically. 
From the Eqs. we get p2] 

0=lH + ^-8a(^O t . (58) 

For the variables <p and £ Eqs. (15511570 constitute a system of second order differential 
equations. It is useful to reduce this system to a first order set of differential equations. 
This can be immediately achieved by introducing new variables, as 

u = <p 2 , v = t (59) 
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Then, we obtain the following first order system 

= -^H 2 + \u + V + 24H 3 v, (60) 

= — (2H + 3H 2 ) + -u-V- 8H 2 v - 16H(H + H 2 )v, (61) 
k 2 

= u + 6Hu + 2V + 2Gv. (62) 

From here it is easy to explore the stability and number of attractors of the physical 
model. To do so we rewrite the system (I60H62I) as 

H = ±GH~ 2 - H 2 , (63) 

ii = —(6Hu + 2V + 2Gv), (64) 

* = 8jp [ i {i k GH ~ 2 - r2 ) + u - ( 16 ™ - 8H >}- ( 65 ) 

Let q = u + 2V. Then this can be put as 

= ~^H 2 + \q + 24H 3 v, (66) 

= —(2H + 3H 2 ) + -q - 2V - 8H 2 v - 16H(H + H 2 )v, (67) 

= q + 6Hq - 12HV + 2Gv, (68) 
and for QS3HB5D 

H = —GH- 2 - H 2 , (69) 

24 K ' 

q = 12HV - GHq - 2Gv, (70) 

* = 8JP [ ¥ { 2h GH ~ 2 ~H 2 ) + q-2V- (16HH - 8H 3 )v]. (71) 
Let us now rewrite Eqs. (I55H57P in terms of the number of e-foldings N. We have 



= —H 3 + \H 2 <? N + V{<j>) + 2AH% N , (72) 
= ^(2HH N + 3H 2 ) + l -E 2 tf N - V{4>) 

-8H% NN - 8H{2HH N + 3H 2 )H£ 1N , (73) 

= + (H N + 3H 2 )<P 2 N + V N + £ 1N G. (74) 

If we introduce new variables 

h = <j) 2 N = H- 2 u, p = £ 1N = H- 1 v, (75) 
the system (I72H74P takes the form 



= -^H 2 + \H 2 h + V{<j>) + 24H*p, (76) 
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= — (2HH N + 3H 2 ) + -H 2 h-V-8H 4 p N -8H(2HH N + 3H 2 )Hp,(77) 
k 2 

= H 2 h N + 2(H N + 3H 2 )h + 2V N + 2Gp. (78) 

To summarize, the original system of equations has been written in several ways that 
can be used to construct solutions. Needless to say, all the corresponding solutions are 
equally valid in every form of the system. Here we will find solutions for the one given 
by Eqs. (I60H62P . To do tat, we rewrite them in terms of the e-folding variable N: 

3 „, 1 
fc2 



= -— H 2 + -u + V + 24:H 3 v, (79) 



= ^(2HH N + 3H 2 ) + ~u-V- 8H 3 v N - 16H 2 (H N + H)v, (80) 
k z 2 

= Hu N + 6Hu + 2HV N + 2Gv. (81) 
This has the following solution 

u = 8H 4 W NN + 8H 3 (3H N - H)W N - —HH^, (82) 

v = HW N} (83) 

V = -4H 4 W NN - {12H 2 H N + 20H 3 )HW N + ^H 2 + ^HH N , (84) 

k k 

where W = W(N) is some function of N. At the same time, 4>,£i are given by 



4> = J ^8H^W NN + 8H(3H N -H)W N - ^(ln H) N dN, (85) 

£ = W. (86) 

Let us now consider several examples. 

(i) As the first one consider the case when W — 0, which corresponds to the Einstein- 
scalar gravity. In this case, we have 

u = -^HH N} v = Q, V = ^H 2 + ^HH N . (87) 

(ii) A second example is the model 

W = u, (88) 
where v = const. Here the solution of the system (4.30)-(4.32) is 

u = -^HH N , v = 0, V = j^H 2 + ^HH N . (89) 
It corresponds to the so-called Einstein-Gauss-Bonnet gravity. 
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(iii) The third example will be 

W = tiN + v (90) 
where /z, v are some constants. Then, the solution of the system (4.30)-(4.32) is 



given by 



2 



u = 8fiH\3H N -H)- —HHn, (91) 

k 

v = ftH(= fxN), (92) 

V = -n(12H 2 H N + 20H 3 )H + ^-H 2 + ^HH N . (93) 

k z k 

The next step is finding the explicit forms of 4>(t), H(t),^i(t), V(t). To this end we 
consider the case when <f> has a kink form, that is, it obeys the (0+l)-dimensional 
Sine-Gordon equation 

= 7 2 sin0, (94) 

where 7= const and 

= Aarctg[e-^- t( % u = 2 7 2 (1 - cos0) = 1 + A j" {t _ to y (95) 
In this case W obeys the equation 

8H 4 W NN + 8H 3 (3H N -H)W N - ^HH N = 2 7 2 (1 - cos0). (96) 
We will solve this equation for the case (4.42) and H = a + (3e mN . There, y satisfies 

ey 2 +5y+a = (12^-8/i)?/ 2 +(12^-16a/i-^r)i/-(M+8/ia 2 ) = 0, (97) 

k l 

which has the solution 

5 ± ^J5 2 - 4ea . , 

y(t) = y{4>) = y t = ^ . (98) 

Hence, we have 

N = N l {t)=N. l {ct>) = -\n y ^ (99) 

and 

H = H(t) = H{(j>) = y/cT+Yi- (100) 
We have thus seen that, from the action (15"5|) . any of the usual cosmologies can be 



achieved. The system of equations (I82ti84p provides a quite simple setup to reproduce 
any kind of cosmological solution, as it has been clearly illustrated with the three 
examples above, where for a given Hubble parameter the scalar Gauss-Bonnet theory is 
constructed. 
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5. Conclusions 

In the present paper, several types of DE cosmologies in modified GB gravity — which can 
be viewed as being inspired by string considerations, Ref.[13] — have been investigated. 
We have studied different kinds of theories, in all of which the GB invariant plays an 
important role in the corresponding equations. First, we have shown that GR plus some 
function of the GB term provides a very powerful theory, where no sort of dark energy 
is actually needed to reproduce the standard ACDM cosmology. 

More general theories have been considered too, as in Sect. Ill, where an action 
depending on a function of the Ricci scalar and GB invariant has been studied. In this 
case the dark energy behavior can also be reproduced, while its extra degrees of freedom 
provide us with a powerful tool to constraint the theory, in order to avoid violation of 
the local gravity tests. 

In the same way, the scalar GB theory studied in the last section adequately 
reproduces any kind of cosmological solution. As a follow up, we have shown, with 
the help of several particular examples corresponding to explicit choices of the functions 
f(G), F(R, G) or the scalar field that, in principle, any cosmic evolution can be obtained 
from these models, what includes the unification of early-time inflation with the late-time 
acceleration coming from astronomical observations. To finish, it has been indicated that 
DE cosmologies in a more general (and complicated) F(R, G) framework — generically 
requiring numerical analysis — can be reconstructed in a similar fashion, too. 

Acknowledgments 

DSG acknowledges a grant from MICINN (Spain), project FIS2006-02842. Part of 
EE's research was performed while on leave at Department of Physics and Astronomy, 
Dartmouth College, 6127 Wilder Laboratory, Hanover, NH 03755, USA. This work was 
also supported by MEC (Spain) grant PR2009-0314, and by AGAUR (Generalitat de 
Catalunya), contract 2009SGR-994. 

References 

[1] Nojiri S., Odintsov S.D., eConf C0602061, 06 (2006) [Int. J. Gcom. Meth. Mod. Phys. 4, 115 

(2007)]. [arXiv:hep-th/0601213] . 
[2] Elizalde E., Nojiri S., Odintsov S.D., Saez-Gomcz D., Faraoni V., Phys. Rev. D 77 106005 (2008) 

[arX"iv:0803.1311j . 

[3] Nojiri S., Odintsov S.D., J. Phys. Conf. Ser. 66 012005 (2007) [arxiv:hep-th/0611071| ; Nojiri S., 
Odintsov S.D., Phys. Rev. D 74 086005 (2006) [arxiv:hep-th7 0608008f~Capozziello S. , Nojiri S., 
Odintsov S.D., Troisi A., Phys. Lett. B 639 135 (2006) [arxiv:astro-ph/0604431| . 

[4] Nojiri S., Odintsov S., Toporensky A., Tretyakov P., Reconstruction and deceleration- acceleration 
transitions in modified gravity |arXiv: 0912.2488 . 

[5] Nojiri S., Odintsov S.D., Saez-Gomez D., Phys. Lett. B 681 74 (2009) |arXiv:0908.1269] . 

[6] Nojiri S., Inverse problem - reconstruction of dark energy models [arXiv:0912.5066 . 

[7] Cognola G., Elizalde E., Odintsov S.D., Tretyakov P., Zerbini S., Phys. Rev. D 79 044001 (2009) 
[arX"iv:0810.4989l. 



ACDM epoch reconstruction from F(R,G) and modified Gauss-Bonnet gravities 



14 



[8] Elizalde E., Saez-Gomez D., Phys. Rev. D 80 044030 (2009) [arXiv:0903.2732 [hep-th]]. 

[9] Nojiri S., Odintsov S.D., Phys. Lett. B 631 1 (2005) | ar xiv:hep-th/0508049| . 

[10] Nojiri S., Odintsov S. D., Sasaki M., Phys. Rev. D 71 123509 (2005) [arXiv:hep-th/0504052l . 
[11] Nojiri S., Odintsov S.D., Gorbunova O.G. J. Phys. A 39 6627 (2006) |arxiv:hep-th/0510183| . 
[12] Cognola G., Elizade E., Nojiri S., Odintsov S.D., Zerbini S., Phys. Rev. D 73 084007 (2006) 
|arxiv:hep-th/0601008| . 

[13] Nojiri S., Odintsov S.D., Sami M., Phys. Rev. D 74 046004 (2006) |arxiv:hep-th/0605039] . 
[14] Zhou S-Y., Copeland E.J., Saffin P.M. JCAP 0907 009 (2009) : arXiv:0903.4610j . 
[15] Goheer N., Goswami R., Peter K. S., Dunsby P., Ananda K., Phys. Rev. D 79 121301 (2009) 
larXiv:0904.2559j . 

[16] Uddin K., Lidsey J.E., Tavakol R. Gen. Rel. Grav. 412725 (2009) |arXiv:0903.0270| . 
[17] Boehmer C.G., Lobo F.S.N., Phys. Rev. D 79 067504 (2009) [arXiv:0902.2982| . 
[18] Alimohammadi M., Ghalee A., Phys. Rev. D 79 063006 (2009) |arXiv:0811.T286] . 
[19] De Felice A., Tsujikawa S., Phys. Rev. D 80 063516 (2009) [arXiv:0907.1830] . 
[20] De Felice A., Tsujikawa S., Phys. Lett. B 675 1 (2009) |arXiv:0810.5712) . 

[21] Cognola G., Elizalde E., Nojiri S., Odintsov S.D., Zerbini S. Phys. Rev. D 75 086002 (2007) 
[heph/0611198| . 

[22] Gurses M., Gauss-Bonnet Gravity with Scalar Field in Four Dimensions [arXiv:0707.0347 . 

[23] Li B., Barrow J.D., Mota D.F., Phys. Rev. D 76 044027 (2007) [arXiv: 0705 .3795] . 

[24] Nojiri S., Odintsov S.D., Trctyakov P.V., Phys. Lett. B 651 224 (2007) larXiv:0704.2520j . 

[25] Bamba K., Odintsov S.D., Sebastiani L., Zerbini S., Finite-time future singularities in modified 

Gauss-Bonnet and .F(i?, G) gravity and singularity avoidance arXiv:091 1.4390 . 
[26] Capozziello S., De Laurentis M., Nojiri S., Odintsov S.D., Phys. Rev. D 79 124007 (2009) 

[arXry:0903.2753| . 

[27] Nojiri S., Odintsov S. D., Phys. Rev. D 68, 123512 (2003) |j hep-th/0307288] . 
[28] Bamba K., Nojiri S., Odintsov S. D., JCAP 0810 045 (2008) |arXiv:0807.2575] 



